The glass transition of glycerol in the volume-temperature plane 
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We assess the relative importance of spatial congestion and lowered temperature in the slowing 
dynamics of supercooled glycerol near the glass transition. We independently vary both volume, V , 
and temperature, T, by applying high pressure and monitor the dynamics by measuring the dielectric 
susceptibility. Our results demonstrate that both variables are control variables of comparable 
importance. However, a generalization of the concept of fragility of a glass-former shows that 
the dynamics are quantitatively more sensitive to fractional changes in V than T. We identify a 
connection between the fragility and a recently proposed density-temperature scaling which indicates 
that this conclusion holds for other liquids and polymers. 
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When a liquid is cooled below its melting point its 
viscosity increases very rapidly with decreasing tempera- 
ture until it turns into a glass, with solid-like properties. 
For instance, glycerol, a typical glass-former, shows a re- 
markable increase in viscosity and relaxation time of 14 
orders of magnitude as it is cooled from 350K to 180K at 
1 atmosphere . There are two broad classes of explana- 
tions for this rapid increase of viscosity. One set of ideas 
identifies density as the crucial variable, arguing that the 
constraints on molecular rearrangements imposed by the 
dense packing in a liquid progressively increase due to 
the thermal contraction that accompanies cooling, until 
finally all motions are arrested in the glassy state. By this 
argument, the driving force for the rapid increase of vis- 
cosity in the case of glycerol would primarily be the 10% 
decrease in volume upon cooling. The second class of ex- 
planations for the slowing dynamics emphasizes the role 
of temperature: lowered temperature renders molecules 
too inactive to move around and surmount the energy 
barriers that impede exploration of their environment. In 
order to resolve the issue of whether it is temperature, T, 
or volume, V , that is the dominant variable in this phe- 
nomenon, experiments are required that independently 
control these variables. However, data of this kind are 
relatively sparse 0. The objective of the experiments 
reported in this Letter is to distinguish unambiguously 
the effects of constrained volume and lowered tempera- 
ture, by using high pressure as a mean of independently 
changing the density of the liquid. 

Glycerol is a widely-studied glass-former with broad 
industrial use. It is also an experimentally convenient 
sample: it has a high dielectric constant and does not 
easily crystallize. Moreover, some cornplementary high 
pressure data exist for glycerol 0; 0> IE la| • These studies 
are principally isothermal experiments in which temper- 
ature is fixed and pressure is varied in discrete steps; as 
such, these data are not optimal for studying tempera- 
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ture dependence. They also appear to disagree with each 
other on the temperature dependence of the relaxation 
frequency Vp near the glass transition temperature, Tg, 
defined as the temperature at which Vp = 0.01 Hz. The 
temperature dependence is often characterized by the 
fragility, m, defined by m = -d{logVp)/d{Tg/T)\T=Tg 
7]. The fragility thus quantifies deviation from the Ar- 
rhenius, or thermally-activated behaviour, and catego- 
rizes the temperature dependence of glass-formers under 
isobaric cooling as strong (small m) and fragile (large 
to). One set of high pressure experiments Q finds that 
the fragility of glycerol increases with high pressure while 
another (5j| finds no pressure dependence of the fragility. 
These results are not necessarily contradictory since they 
are taken in different frequency regimes, involving differ- 
ent degrees of extrapolation to infer the behaviour at 
Tg. The present measurements go down to 0.01 Hz, ob- 
viating the need to extrapolate down to Tg in order to 
infer fragility. We have also performed an isobaric ex- 
periment by controlling pressure, P, and studying the 
temperature dependence in detail. We have restricted 
ourselves to pressures below 1 GPa and taken data at 
several closely separated pressures to study changes in 
dynamics as functions of P and V as well. 

We employ a standard Teflon cap technique, widely 
used in high-pressure experiments The pressure is 
inferred from the resistivity of a manganin coil immersed 
in the liquid. Manganin is an alloy whose resistivity as 
a function of P and T is known J^J . After the liquid is 
squeezed to a high pressure at room temperature, we cool 
it down to low temperatures. The dynamics of the liq- 
uid is studied by measuring the capacitance of a coaxial 
capacitor immersed in the liquid. We obtain the com- 
plex dielectric susceptibility e(j/) — e'iv) -I- i(" [y) for v 
ranging from 10~^ to 10^ Hz at each temperature. An 
example of the imaginary part of the response, (" {v) is 
displayed in the inset of Fig. ^ We extract the charac- 
teristic relaxation frequency, j/p, defined as the frequency 
of the peak in e"{iy), at each P and T. In addition to the 
characteristic frequency, e^v) also yields the distribution 
of relaxation times and the ionic conductivity [3] . In this 
Letter we focus on the dependence of Vp on T, V, and P. 
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FIG. 1: Inset: e" vs. v at 217 K and 340 MPa. The line 
is Davidson-Cole fit: e = too + (eo — eoo)/(l — iv/faY . The 
relaxation frequency, Vp is defined to be the peak frequency. 
Main figure: Isobaric plot of Vp vs. T. The circles repre- 
sent data at 83, 164, 212, 278, 358, 392, 527, 539, 560, 689, 
876 ± 4 MPa. The crosses are 1 atmosphere data l]. The 
solid fines are Cohen-Grest fits: = A^T -b/[T - Tc + 

\/ {T ~ TcY + CT. A, b, and C are pressure-independent. 
Tc — VaP where Va is a volume parameter. The theory pre- 
dicts that the molecular volume, Vm ~ 4.605(6/C) Va whence 
Vm = 1700A^ 

In Fig. n we display i/p vs. T for 12 different pressures. 
As is evident from the figure, pressure has a strong ef- 
fect on relaxation dynamics: our maximum pressure of 
about 1 GPa decreases the relaxation frequency at a fixed 
temperature by about seven decades. 

TurnbuU and Cohen 11] proposed a theory in which 
the slowing dynamics in a supercooled liquid is due to 
reduced free volume at low temperature. Free-volume 
theory predicts that the viscosity at a fixed temperature 
is given by r/ = exp(l^/V/) where Vm is the molecular 
volume and Vf the free volume. The temperature depen- 
dence of rj was computed by Cohen and Crest who 
derived a formula for the viscosity at a fixed pressure. 
Using the Einstein-Debye relation, i/p — BT/i], where B 
is a constant with the dimension of volume, we compare 
this prediction to the data in Fig. 2] While the fit is good, 
the fit parameters lead to a prediction for a mean molec- 
ular volume of ITOOA^, in excess of length scales from 
neutron scattering J^] or even of a simple estimation of 
123 A^ from the mass density and molar mass. 

A much more direct way to test free-volume theory 
is to study the dynamics when the volume - and there- 
fore the free volume - is held fixed. Usingan empirical 
equation of state for glycerol we calculate 'q] the volume 
corresponding to each T and P in Fig. ^ This yields 
the isochoric and isothermal plots shown in Fig.|21 Free- 
volume theory predicts [Tlj that the viscosity varies only 
weakly with temperature {rj ^ VT) at constant volume. 



Fig. 121 a) shows that the isochoric change of Vp with re- 
spect to temperature is a factor of 10^, much faster than 
the factor of 1.1 predicted by the free volume model. 

It is clear from Fig. I2a) that even when the volume 
is constant, the dynamics remains strongly temperature 
dependent. Indeed, we observe from Fig. ^a) that the 
1 atmosphere isobaric curve has only a slightly stronger 
temperature dependence compared to nearby isochoric 
curves. This means that there is only a small volume 
contribution to the slowing dynamics when the liquid is 
supercooled along the 1 atmosphere isobar. This obser- 
vation has led Ferrer et al. to conclude that temper- 
ature is overwhelmingly more important than volume. 
Does this more generally imply that temperature is the 
dominant variable in the glass transition and that the 
effects of volume are secondary? 

To address this question we isolate the role of volume 
that of temperature by in an isothermal plot, such as the 
one shown in Fig.|2Ib) in which volume is varied, holding 
T fixed. We observe in Fig. Efb) that I'p changes strongly 
with l/V, faster than exp{l/V). There appears to be 
greater variation in Vp in the isochoric cases, however 
we point out that l/V increases only 14% in Fig. |2Ib) 
whereas 1/T increases by 54% in Fig. [2Ia). 

To assess the relative contribution of V and T it 
is therefore important to take into account equiva- 
lent fractional changes in these variables. To this 
end we compare A(\ogVp)\v for some AT/T with 
A(\ogi'p)\T for the same AT/T. This leads us to con- 
sider the variables, Mt = —T{d/dT) \ogVp\v and My = 
—V{d/dV)\ogVp\T- When evaluated at Vp = 0.01 Hz, 
Mv is the isothermal fragility with respect to V and 
Mt the isochoric fragility with respect to T. We call 
these quantities the generalized fragilities by analogy to 
the traditional fragility, to, which compares glass formers 
with different Tg because to measures the rate of slowing 
dynamics only after scaling out the material-dependent 
Tg. Likewise My and Mt compares the rate of change 
of Vp in V and T directions only after scaling out the 
relevant variables. Furthermore, the variables My and 
Mt are manifestly symmetric in V and T . 

We display in Fig. Ola) the vector field (Mt, My) plot- 
ted in the V-T plane. At all places where we can evaluate 
Mt and My, they are comparable: 1.2 < {My/MT) < 
1.7. Thus the effects on relaxation frequency of equal 
fractional changes in volume or temperature are com- 
parable over the entire experimental range, and neither 
variable can be neglected. While neither variable is over- 
whelmingly dominant, it is clear that My > Mt over the 
experimentally accessible region. Thus relaxation fre- 
quencies are more sensitive to fractional changes in V 
than T everywhere in the V-T plane. In particular, the 
data indicate that this is also true for the 1 atmosphere 
isobaric glass transition at Tg — 192 K, Vg = 0.959. 
While we can only directly measure Mr at (192,0.959) a 
linear extrapolation of My to the smallest temperature 
in Fig.l^a) shows that My/MT = 1-58. 

There are two opposing factors affecting the fragilities 
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FIG. 2: (a) Isochoric plot. The solid lines are Vogel-Fulcher-Tammann-Hesse (VFTH) fits: Vp = 

I'oo exp[i?/(ro — T)]. The heavy line shows 1 atmosphere isobaric data = 1 at 273.15 K, 1 atmosphere. 

Data are shown for y=0. 970,0. 951,0. 942, 0.931,0. 921,0. 916,0. 900,0. 895,0. 895,0.878, 0.860 ± 0.3%. (b) Isothermal 
plot. These data are obtained from isothermal cuts to Fig. The temperatures from left to right are: 
214,216,218,220,222,224,226,228,230,232,234,236,238,240,242 K. The solid lines are Doolittle fits: Vp = v^c exp[D/(l/o - V)]. 




FIG. 3: (a) Field of fragility vector {Mt,Mv)- The scale bar shows the size of the vectors. The dashed line is a 45° reference 
line showing that Mv > Mt for all vectors. We can evaluate only one component for the grey vectors. The solid curves from 
top left to bottom right show the glass transition line of constant relaxation frequency 0.01 Hz; the horizontal and vertical 
grey vectors at this glass transition line are the isochoric and isothermal fragilities respectively, (b) Test of V-T scaling, 
\ogVp = /{VT). The best possible collapse of data is achieved with 7 = 1.4 



as one moves along the glass line towards high densities: 
increasing T reduces fragility while decreasing V makes 
it bigger. It is clear from Fig. |3Ia) that the latter effect 
dominates, inasmuch isochoric and isothermal fragility 
both increase along the glass transition line towards lower 
T or lower V. The increase of isochoric fragility, Mr, 
with density is an additional indicator of greater sensi- 
tivity to changes in V than T. We have observed that 
the width of the relaxation too, increases more rapidly in 



the V direction than the T direction @. 

Are these results specific to glycerol? We address 
this issue by means of a density-temperature scaling 
suggested by Casalini and Roland and also Alba- 
Simionesco et al. |0. For several liquids and polymers 
they found that logz/p = f{V^T) where 7 and the func- 
tion / are material dependent. Next, we observe that 
this scaling exponent is equal to the ratio of the general- 
ized fragilities: 7 = My /Mt- This scaling for our data is 
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FIG. 4: Glass transition in VT plane. The solid lines are 
constant frequency contours. The error bars are of compara- 
ble size to the symbols. The solid circles are To from isochoric 
VFTH fits and the open circles Vb from isothermal fits to the 
Doolittle equation to data sets that extend to 0.01 Hz. The 
open squares are taken from Ref Q. The dashed lines involve 
greater extrapolations of the data. 

shown in the Fig.|3tb) with 7 = 1.4; this is less than 1.8 
obtained usine; the data of Cook et al. by Ref. [13 ■ 
The scaling in Fig.|3fb) is by no means perfect, reflecting 
the fact that Mv/Mt is not constant. For a power-law 
potential, V{r) ~ r~^^, the scaling is exact. However, in 
the more complicated situation of a molecular liquid, one 
should not be surprised if the scaling is inexact because 
changes in density - and consequently the intermolecular 
spacing - will emphasize different parts of the potential. 
Nevertheless, to the extent that this scaling is even ap- 
proximately successful over a limited range of data, 7 > 1 
implies that My > Mt- There are 25 liquids with 7 > 1 
and only one (sorbitol) with 7 < 1 0] ■ Therefore 7 > 1 is 
observed for a wide range of liquids. The greater sensitiv- 
ity of relaxation dynamics to fractional volume changes 



is thus not limited to glycerol, and appears to be typical 
rather than exceptional. 

We summarize the new data for glycerol in Fig. ^ as 
a set of constant relaxation frequency curves in the V- 
T plane. The glass transition line {vp =0.01 Hz) is 
formed by Tg and Vg extracted from Fig. [3 The con- 
stant frequency curves get denser towards lower frequen- 
cies. Apart from the constant frequency lines in the ex- 
perimentally accessible regime we also indicate in Fig. 0] 
extrapolations to the zero frequency limit. The solid cir- 
cles represent Tq from VFTH fits to isochoric data and 
the open circles Vb from Doolittle fits to isothermal data. 
While the existence of a finite temperature or finite vol- 
ume structural arrest depend solely on extrapolations, 
it is worth pointing out that these two independent ex- 
trapolations support each other in that they appear to 
produce a continuous (Vb,Tb) curve. With the available 
data the value of Tq appears to be increasing in an un- 
bounded way at high densities. This agrees with the ex- 
pectation that an extremely compressed system will re- 
semble a hard sphere liquid, where Tq ^ 00. This limit 
is approached along a concave curve, as in attractive col- 
loidal systems 18]. Clearly, even broader explorations in 
the V-T plane of both dynamics and of the Kauzmann 
entropy crisis would be useful in assessing the role of Tq 
and Vq. We end by recapitulating our principal result: 
neither axis on this plane may be neglected in under- 
standing supercooled glycerol. 
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